Checking the Winter timeseries for when the difference between the top and the bottom is greater than -0.1. The red line “inverse” is the 1m thermistor minus the 9m thermistor while “inverse 2” is the surface thermistor minus the 9m thermistor. The red dotted line is at -0.1 because the accuracy of the thermistors is 0.075.
After looking at the the timeseries let’s break it down for each period and look more closely at the temperature profiles, disolved oxygen profiles, and images from the camera at the catwalk.
profile_plot(DATA=FCR_Catwalk_EDI_2018_2023,
start_time = "2022-12-22 00:00:00",
end_time = "2023-01-05 00:00:00",
Varibles = starts_with("Thermis"))
## Warning: Using an external vector in selections was deprecated in tidyselect 1.1.0.
## ℹ Please use `all_of()` or `any_of()` instead.
## # Was:
## data %>% select(Varibles)
##
## # Now:
## data %>% select(all_of(Varibles))
##
## See <https://tidyselect.r-lib.org/reference/faq-external-vector.html>.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
# DO
profile_plot(DATA=FCR_Catwalk_EDI_2018_2023,
start_time = "2022-12-22 00:00:00",
end_time = "2023-01-05 00:00:00",
Varibles = c("RDOsat_percent_5", "RDOsat_percent_9", "EXODOsat_percent_1"))
Pretty clear inverse stratification from about Dec. 23 to January 2nd.
It looks like we see ice on the camera on Dec. 25th
Thin layer of ice on the reservoir
profile_plot(DATA=FCR_Catwalk_EDI_2018_2023,
start_time = "2022-12-31 00:00:00",
end_time = "2023-01-05 00:00:00",
Varibles = starts_with("Thermis"))
# DO
profile_plot(DATA=FCR_Catwalk_EDI_2018_2023,
start_time = "2022-12-22 00:00:00",
end_time = "2023-01-05 00:00:00",
Varibles = c("RDOsat_percent_5", "RDOsat_percent_9", "EXODOsat_percent_1"))
It looks like there is no longer inverse stratification on some time on January 2nd. Let’s look at the pictures
There could be faint ice in the backaground
Don’t see any ice in this picture and getting warm
profile_plot(DATA=FCR_Catwalk_EDI_2018_2023,
start_time = "2023-01-15 00:00:00",
end_time = "2023-01-18 00:00:00",
Varibles = starts_with("Thermis"))
# DO
profile_plot(DATA=FCR_Catwalk_EDI_2018_2023,
start_time = "2023-01-15 00:00:00",
end_time = "2023-01-18 00:00:00",
Varibles = c("RDOsat_percent_5", "RDOsat_percent_9", "EXODOsat_percent_1"))
Could be ice near the far shore
Don’t see any ice in this picture and getting warm
Raining and no ice